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A  review  of  the  research  chat  has  been  carried  out  under  the  sponsorship  of  Air 
Force  Office  of  Scientific  Research  grant  number  AFOSR-90-01 79  is  given.  Details 
are  given  about  two  research  programs,  namely  flow  control  and  superconductivity. 

In  both  instances  we  have  modeled  realistic  physical  problems,  we  have  performed 
rigorous  analysis  on  the  mathematical  models,  we  have  developed  approximation 
schemes  which  we  have  al.so  analyzed,  and  we  have  written  computer  codes 
Implementing  our  computational  method.  In  both  Instances,  were  the  first  to  carry 
out  such  a  comprehensive  research  program,  and  in  so  doing,  have  developed  robust 
and  useful  algorithms  that  also  embarked  on  a  joint  research  program  with  personnel 
at  AEDC  facility  in  Tullahoma.  Tennessee.  We  also  briefly  discuss  other 
accomplishments  related  to  this  grant.  Including  other  research  project  invited 
talks,  journal  articles,  etc. 
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ABSTRACT 


A  review  of  the  research  that  has  been  carried  out  under  the  sponsorship  of  Air  Force  Office  of 
Scientific  Research  grant  number  AFOSR-9O-0179  is  given.  Details  are  given  about  two  research 
programs,  namely  flow  control  and  superconductivity.  In  both  instances  we  have  modeled  realistic 
physical  problems,  we  have  performed  rigorous  analyses  on  the  mathematical  models,  we  have  de¬ 
veloped  approximation  schemes  which  we  have  also  analyzed,  and  we  have  written  computer  codes 
implementing  our  computational  method.  In  both  instances,  we  were  the  first  to  carry  out  such  a 
comprehensive  research  program,  and  in  so  doing,  have  developed  robust  and  useful  algorithms  that 
may  be  applied  to  practical  problems.  In  connection  with  our  research  in  flow  optimization,  we  have 
also  embarked  on  a  joint  research  program  with  personnel  at  the  .\£DC  facility  in  TuUahoma,  Ten¬ 
nessee.  We  also  briefly  discuss  other  accomplishments  related  to  this  grant,  including  other  research 
projects,  invited  talks,  journal  articles,  etc. 
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I  -  REVIEW  OF  RESEARCH  ACCOMPLISHMENTS 
UNDER  GRANT  NUMBER  AFOSR- 90-01 79 


We  begin  by  reviewing  some  of  the  research  that  has  been  completed  and  that  has  been  supported 
by  AFOSR  under  Grant  Number  AFOSri-90-OlT9  For  the  sake  of  brevity,  we  wili  not  go  into  great 
detail  in  the  following  diacuasion;  further  information  concerning  these  topics  can  be  gauned  bom 
the  appropriate  references  listed  m  §III.  .Also,  again  for  the  sake  of  brevity,  will  only  diacuss  two 
of  the  areas  of  research  with  which  we  have  been  concerned  Other  research  projects  are  bneUy 
discussed  in  §11. 


1.  Numerical  Solution  of  Flow  Control  Problems 

Over  the  past  few  years  we  have  undertaken  a  comprehensis'e  research  program  addreaating  vanous 
issues  connected  with  3ow  control  and  optimization  problems  Our  efforts  in  this  direction  have 
focused  on 

•  building  mathematical  models  of  the  physical  problems,  invoking  a  mimmum  of  assump¬ 
tions  about  the  physical  phenomena; 

•  rigorously  analyzing  the  mathematical  models,  for  example,  to  study  the  existence  and 
regularity  of  solutions,  to  verify  the  existence  of  Lagrange  multipliers  to  enforce  ojnstraints, 
and  to  derive  necessary  conditions  that  optimal  controls  must  satisfy; 

•  constructing  and  analyzing  discretization  methods  for  determining  approximate  solutions 
of  the  optimal  control  and  optimization  problems,  including  a  rigorous  derivation  of  error 
estimates;  and 

•  developing  computer  codes  implementing  our  discretizauion  algorithms,  first  for  the  purpose 
of  showing  the  efficacy  of  these  methods,  and  ultimately,  to  solve  problems  of  practical 
interest. 

Our  success  in  carrying  out  the  program  is  evidenced  in  various  ways.  V/c  have  published,  or  have 
had  accepted,  numerous  articles  in  top-of-the-line  journals;  we  have  been  asked  to  write  chapters  for 
books;  we  have  been  invited  to  lecture  at  many  international  conferences  both  of  the  mathematical 
aund  engineering  persuasions;  we  were  asked  to  organize  and  give  the  main  talk  at  a  workshop  on  Sow 
control  held  in  November,  1992  at  the  Institute  for  Mathematics  and  its  Applications  in  Minnesota; 
and,  perhaps  most  important,  our  ideas  and  algorithms  are  been  used  and  implemented  by  engineers 
interested  in  flow  control  and  optimization  problems.  Details  concerning  the  papers,  invited  talks, 
and  the  Minnesota  workshop  may  be  found  in  §III;  a  discussion  of  a  setting  wherein  our  ideas  are 
being  implemented  is  discussed  in  §1.1. 

Our  work  on  flow  control  deals  with  two  different  classes  of  methods  of  solution.  In  §1.1  we 
discuss  our  work  on  methods  that  use  sensitivities  within  an  optimization  algorithm,  while  in  §1.2 
we  deal  with  the  adjoint,  or  co-state,  or  Lagrange  multiplier  approach. 

1.1.  Design  of  the  AEDC  free-jet  test  facility 

An  example  of  our  efforts  in  using  sensitivities  in  connection  with  flow  control  problems  is  the 
work  we  are  doing  in  collaboration  with  engineers  at  the  AEDC  facility  in  Tillahoma,  Tennessee. 
Our  collaborators  in  this  project  include  John  Burns,  Gene  Cliff,  and  Jeff  Borggard  of  Virginia  Tech. 
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aad  Karl  Kaeile  aad  Doaaid  Todd  of  A£DC 

One  of  the  primary  missions  of  the  AEDC  facility  is  to  oast  aixcralt  engines.  Indeed,  within 
the  United  States,  they  have  the  unique  capability  of  testing  fuii-s«e  enginw  within  the  twitraded 
environment  of  a  wind  tunnel.  Ideally,  one  would  like  to  have  the  Sow  entering  the  engine  m  the 
test  facility  to  be  exactly  the  same  as  that  entering  the  engine  if  it  were  attached  to  an  aircraft  in 
Sight.  Thus,  referring  to  the  sketch  in  the  Sgure,  we  would  like  the  Sow  at  a  reference  plane  aliead 
of  the  engine  in  the  test  facility  to  be  the  same  as  that  in  a  similarly  situated  plane  ahead  of  tlie 
engine  in  Sight.  In  order  to  match  the  Sight  conditions  at  the  reference  plane  (the  test  fsciiily  is 
not  large  enough  to  accommodate  the  whole  aircraft)  one  is  free  to  choose  the  forebody  shape  and 
certain  inSow  parameters  such  as  the  total  pressure  and  .Mach  number:  see  the  Sgure. 


REFERENCE 

PLANE 


TEST  CELL  WALL 


Thus  the  optimization  problem  is  easily  stated:  determine  the  inSow  parameters  and  the  fore¬ 
body  shape  that  minimizes  the  difference  between  the  Sow  at  the  reference  plane  in  the  test  facility 
and  that  in  Sight. 

In  order  to  explain  our  ideas  in  connection  with  this  problem,  we  introduce  some  notation.  The 
mathematical  structure  of  our  problem  is  similar  to  that  of  any  optimization  or  control  problem. 
First  we  have  certain  state  variables,  i.e.,  variables  that  serve  to  describe  the  state  of  the  system.  In 
our  context  these  are  given  by  the  prewure  p,  the  internal  energy  c,  the  density  p,  and  the  momentum 
Seld  pu,  where  u  denotes  the  velocity.  We  collectively  denote  the  state  variables  by  <i>.  Next,  we 
have  control  variables  or  design  parameters  which  we  collectively  denote  by  g.  In  our  context  these 
include  certain  inflow  conditions  such  as  the  total  pressure  and  Mach  number,  as  well  as  a  finite  set 
of  parameters  that  determine  the  shape  of  the  forebody.  (Typically,  the  forebody  is  defined  as  a 
combination  of  Bezier  curves  that  are  determined  by  a  finite  set  of  parameters.)  The  next  ingredient 
is  an  objective  or  cost  functioned  J  that  in  general  depends  on  the  state  and  may  also  depend  on 
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the  oontrols.  In  our  case  we  c«m  uae 


J{xx)^\[  (1) 

^  Jrnf- 

where  Vnp  denotes  the  reference  plane  and  U/  denotes  the  velocity  Eeid  at  the  reference  plane 
in  flight.  (The  latter  is  obtained  through  experimental  tneaaurementa.)  The  final  ingredients  are 
constnunts  which  in  our  context  we  denote  by  F(ip,  3)  =  0.  These  are  merely  the  flow  equations,  or 
in  a  practical  setting,  discretized  versions  of  the  flow  equations.  Then,  the  problem  in  hand  m  the 
following  mimmuation  problem:  find  controls  g  and  states  <9  such  that  Jiu)  is  minimized,  subject 
to  F{<l>,g}  =  0. 

The  approach  used  at  A£DC  for  solving  this  optimization  problem  relies  on  the  Gauss-Sewlon 
method; 

Start  with  an  initial  guesa  go  for  the  design  parameters. 

For  Ic  =  0, 1,2, . . . , 

1.  use  the  discretized  Sow  equations  to  And  a  state  Ok  such  that  F(Ok.Sk}  =  0; 

2.  compute  the  sensitivities  ido/dg)jg^ 

3.  sol-ve  for  the  increment  5  from 


A 


/  do  do 

\dg  dg 


\ 

dg  \,J  ' 


where  A(-,  )  and  B{-,  )  are  appropriately  defined  bilinear  forms  and  r  denotes  an  appro¬ 
priate  residual. 

4.  Set  Qk+i  =  +  <5 

The  iteration  is  terminated  whenever  or  perhaps  iJ'fu*)  -  J'(u*_i)|,  is  smaller  than  some 
prescribed  tolerance. 

Two  key  observations  concerning  the  above  algorithm  are  that  Step  1 ,  the  state  calculation,  has 
been  carried  out  using  PARC-codes  developed  at  AEDC  and  that  for  Step  2.  the  calculation  of  the 
sensitivities,  finite  difference  approximations 


d<t>\  _«>(pt)-«{5) 

■5-  ~ - :: —  V*; 

09  9k -9 

are  used,  where  g  is  some  value  close  to  gk,  and  <^{g)  satisfies  the  constraint  equations  F{0,9)  -  0. 
Thus,  using  (2),  Steps  1  and  2  require  multiple  state  calculations;  indeed,  one  needs  iV  4-  1  state 
calculations,  where  N  denotes  the  total  number  of  design  parameters.  The  necessity  of  multiple 
state  computations  make  this  approach  prohibitively  expensive,  especially  in  three-dimensions. 

Our  first  task  was  to  make  the  calculation  of  Step  2  more  eflflcient.  We  were  constrained, 
for  practical  purposes,  to  use  existing  PARC-codes,  or  small  variations  of  th^.  In  this  context, 
we  suggested  that  the  sensitivities  should  be  computed  directly  as  solutions  of  the  state  equations 
Fi4»,g)  ~  0  differentiated  with  respect  to  the  design  parameter,  i.e., 


(—1  ^ 

do  _ 

dF 

dg  g^~ 

dg 

(3) 


Note  that  and  (dF depend  only  on  gk  and  <f>k,  so  that  they  may  be 

evaluated  after  Step  1.  Then,  one  may  solve  the  linear  system  (3)  for  the  sensitivities  {dif)/dg)\g^. 
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We  have  determined  how  to  compute  the  sensitivities  within  the  PARC  environment,  and  v«tri- 
anta  of  a  2-D  PARC  code  have  been  developed  both  at  Virginia  Tech  and  AEDC  which  in  fact 
compute  the  sensitivities.  Currently  these  codes  have  been  run  for  the  inflow  and  forebody  design 
parameters.  Results  are  very  encouraging,  e  g.,  savings  of  factors  of  5  to  10  over  the  use  of  the 
difference  approximation  (2)  have  been  demonstrated. 

1.2.  Anaiysij  and  approximatton  of  adjoint  methods  for  flow  control 

.As  was  mentioned  above,  the  structui..  of  a  flow  control  or  optimization  problem  involves 
objectives  chat  express  why  one  wants  to  control  the  flow.  Such  goals  are  usually  expressed  in  terms 
of  the  minimization  of  a  cost  functional  Also,  we  have  constraints  that  determine  what  kind  of  Sow 
one  is  interested  in.  The  nature  of  the  flow  is  expressed  in  terms  of  a  specific  set  of  flow  equations, 
e.g.,  Euler,  Navier-Stokes,  stationary,  time  dependent,  incompressible,  compressible,  etc.  Finally, 
one  has  available  control  mechanisms  or  design  parameters  that  are  to  be  used  to  meet  the  objective. 
Controls  are  expressed  in  terms  of  unknown  data  in  the  problem  specification. 

In  our  own  work  on  adjoint  methods,  we  have  considered  a  targe  variety  of  objectives  and 
control  mechanisms  in  Che  context  of  incompressible  viscous  Sows.  We  have  rigorously  analyzed 
these  methods  with  regard  to  the  existence  and  regularity  of  optimal  solutions  and  the  existence 
of  Lagramge  multipliers,  we  have  derived  optimality  systems  that  determine  optimal  cont.’^ols  .and 
states,  we  have  defined  and  analyzed  finite  element  algorithms,  and  we  have  developed  codes  that 
implement  the  algorithms. 

Among  the  objectives  that  we  have  treated  are  the  following;  this  list  is  not  exhaustive,  but  is 
merely  representative.  First,  we  have  flow  tracking  wherein  one  wisher  the  velocity  field  to  be  as 
close  as  possible  to  a  prescribed  field.  If  u  and  Uj  denote  the  velocity  field  and  a  prescribed  velocity 
field,  respectively,  then  we  want  to  control  the  flow  so  that  u  is  “close”  to  U.^.  For  example,  one 
can  minimize  the  functional 

J,(u)  =  i  f  \u-Vs\*dn 
Jn 

where  fl  denotes  the  flow  domain.  A  second  objective  we  have  treated  is  viscous  drag  minimization. 
This  can  be  accomplished  by  minimizing  the  integral  of  the  dissipation  function,  i.e., 

J2(u)  =  ^^|(Vu)  +  (Vu)‘i’<fQ. 

where  (i  denotes  the  viscosity  coefficient.  Another  important  objective  is  the  avoidance  of  hot  spots, 
i.e.,  places  where  temperature  peaks  occur,  along  bounding  surfaces  since  often  such  phenomena  lead 
to  meltdown  or  to  flexural  failures.  Such  difficulties  may  be  avoided  by  minimizing  the  functional 

MT)=  f  ir-r.i^dT, 

Jr. 

where  T  denotes  the  temperature,  F,  the  portion  of  the  boundary  along  which  one  would  like  to 
avoid  hot  spots,  and  T,  a  desired  temperature  distribution. 

One  of  the  control  mechanisms  that  we  have  considered  is  the  velocity,  or  mass  flux,  along 
portions  of  the  boundary,  i.e.,  injection  or  suction  of  fluid.  Thus,  if  Fe  denotes  the  portion  of  the 
boundary  covered  by  orifices,  we  would  seek  a  control  g  such  that  a  desired  objective  functional  is 
minimized,  subject  to  the  appropriate  flow  equations,  and  also  such  that 

u  =  g  on  Fc . 
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A  second  control  mechaxusm  chat  we  have  considered  is  heating  or  caoiing  along  btoundmg  surface  $. 
For  example,  one  could  seek  a  control  q  such  that  a  desired  objective  functional  is  minitruxed,  subjec*? 
to  the  appropriate  how  equations,  and  also  such  that 


dn 


=  <? 


on  fc . 


where  denotes  the  portion  of  the  boundary  along  which  one  allows  the  control  to  act  and  O/iM 
denotes  the  normal  derivative  at  the  boundary.  A  third  type  of  control  is  a  dtstnbuted  i^antwL 
e.g.,  one  could  try  to  effect  control  through  the  body  force  in  the  .Vavier-Stokes  equation  or  a  heat 
source  in  the  energy  equation.  Thus,  one  would  seek  a  control,  defined  on  the  fiow  domain  Q  or 
on  a  portion  of  n,  such  that  some  functional  is  minimized  and  subject  to  the  appropriate  flow 
equations.  Physically,  one  may  effect  such  control  through  a  magnetic  field  acting  on  an  ionized 
fluid  or  electrically  conducting  fluid  in  the  first  case,  or  through  radiation  mechanisms  or  a  targeted 
laser  beam  in  the  latter  case.  Recently  we  have  begun  studying,  in  the  context  of  adjoint  methods. 
shape  controls,  in  this  case  control  is  effected  by  adjusting  the  shape  of  the  flow  domain.  The  shape 
of  the  flow  domain  may  be  changed  in  many  ways.  For  example,  one  could  iLse  leading  .and/or 
trailing  edge  flaps,  or  movable  walls,  or  rudders,  or  propeller  pitch.  A  related  problem  is  the  optimal 
design  problem.  Here,  we  want  to  choose  a  flow  domain,  e.g.,  the  exterior  of  an  airfoil,  so  that  some 
objective  is  achieved. 

1.3.  Computational  results  for  a  cooling  control  problem 

In  order  to  illustrate  the  type  of  results  that  we  can  obtain  using  adjoint  methods,  we  consider 
the  problem  of  avoiding  hot  spots  along  the  top  wall  of  a  rectangular  channel.  We  assume  that 
the  flow  is  incompressible  and  convection  driven  so  that  buoyancy  effects  can  be  neglected,  and 
thus  temperature  effects  on  the  mechanical  properties  of  the  Bow,  i.e.,  the  velocity  and  pressure, 
are  negligible.  We  are  interested  in  the  design  of  heating  and  cooling  controls  such  that  hot  spots 
are  avoided  along  a  portion  of  the  boundary,  and  thus  we  assume  that  the  flow  is  stationary.  As 
a  result  of  our  assumptions  about  the  flow,  the  state  variables,  i.e.,  the  velocity  u.  pressure  p,  and 
temperature  T,  are  required  to  satisfy  the  continuity  equation,  the  Navier-Stokes  equations,  and  the 
energy  equation. 

In  order  to  avoid  hot  spots,  we  minimize  the  functional 

=  ir-T.i"dr  +  ^  \g\^dr,  (4) 

where  T  is  required  to  satisfy  the  flow  equations  along  with  Naundary  conditions.  In  particular, 
along  Fc,  the  part  of  the  boundary  along  which  controls  are  allowed  to  act,  we  have  that  dTfdn  ~  g. 
r,  denotes  the  part  of  the  boundary  along  which  we  wish  to  match  the  temperature  T  to  &  desired 
temperature  distribution  T,.  The  above  functional  has  been  penalized  by  a  norm  of  the  control; 
this  is  necessary  since  we  are  not  placing  any  a  priori  constraints  on  the  size  of  the  control.  The 
pa,  jjneter  6  is  used  to  adjust  the  relative  sizes  of  the  two  terms  contributing  to  the  functional. 

The  flow  domain  is  the  rectangular  channel  depicted  in  the  accompanying  figure.  All  computa¬ 
tions  were  carried  out  using  piecewise  linear  finite  elements  on  a  triangular  mesh. 
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la  che  Qext  igure  we  show  'he  uncoacroUed  cemperature  ;n  che  interior  aa  'weil  as  along  the  top 
(wadi)  and  right  outflow)  boundanea.  One  sees  that  there  are  many  peaks  along  the  top  boundary, 
amd  shat  there  are  peaks  in  the  temperature  at  the  top  and  bottom  of  the  outhow  ’boundary. 


1 


Top:  level  lines  of  the  uncontrolled  temperature; 

Bottom-left:  uncontrolled  temperature  along  the  top  wall  Ps; 
Bottom-right;  uncontro.  ^  temperature  along  the 
right  (outflow)  boundary  T.vi. 
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Next,  we  have  a  heat  flux  control  acting  adong  r^,  the  cop  boundary,  so  that 


dn 


-'J 


on  , 


where  g  denotes  the  control.  We  also  have  chat  F,  =  Tc,  the  top  boundary,  with  the  target  temper¬ 
ature  along  r,  appearing  in  the  functional  (4)  taken  to  be 


T,  =  cain{3.5,  li.STSx  -^1}  on  F,  =  . 


(5) 


This  function  ensures  a  continuous  transition  with  respect  to  the  indow  temperature  which  h  set  to 
T  =  1.  In  the  next  figure  we  display  the  optimal  temperature  T  and  co-state  in  the  interior  and  the 
optimal  temperature  T  and  control  g  along  the  top  wail  F,  =  F-.  We  see  from  the  graph  of  T  aiong 
r,  that  we  have  very  effectively  matched  the  gi'^Q  distribution  T,  given  by  (S).  The  control  in  this 
case  involves  both  cooling  and  heating. 


'  -V  •w 


Control  and  target  temperature  on  the  top  wall  . 

Top;  level  lines  of  the  optimal  temperature; 

Left-center;  optimal  temperature  along  the  top  wail  F,  =  Fci 
Right-center;  optimal  control  along  the  top  wall  Fc: 

Bottom;  level  lines  of  the  optimal  co-state. 
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Next,  we  choose  F,  =  F.vi,  the  right  or  outflow  boundary,  so  that  F,  F^;  we  aow  have  the 
target  temperature 

Ts  =  3.5  on  F,  =  F.vi .  (6) 

.A.gain,  in  the  next  Sgure,  we  display  the  optimal  temperature  T  and  co-state  in  the  iacenor.  the 
optimal  temperature  T  along  the  outflow  'ooundary  Fj  =  F.vi,  and  control  g  along  the  top  wail  Fc- 
Once  again,  an  excellent  job  is  done  of  matching  to  the  prescribed  temperature  distribution  ;6). 
Note  chat  the  optimal  control  requires  cooling  along  the  cop  boundary  only  m  the  vicuur/  of  the 
outflow  boundeuy,  i.e.,  near  the  location  at  which  we  are  crying  to  match  the  flow  temperature  to 
that  given  by  (6). 


:3  ■ 


!  1 


-  ;.;5 


Control  on  the  top  wall  Fc  and  target  temperature  on  the  right  (outflow) 
boundary  F,  =  F.vi- 

Top:  level  lines  of  the  optimal  temperature; 

Left-center:  optimal  temperature  along  the  outflow  boundary  F,  =  F.vi; 
Right-cemsr;  optimal  control  along  the  top  wall  F^; 

Bottom:  level  lines  of  the  optimal  co-state. 
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2.  Analysis  and  Approximation  of  Macroscopic  Models  for  Superconductivity 

Superconductivity  is  one  of  the  Grand  Challenges  identified  as  being  crucial  to  future  economic 
prosperity  and  scientific  leadership.  Much  of  the  research  in  superconductivity  addresses  the  mi¬ 
crostructure  of  superconductors;  one  goal  of  that  research  is  to  predict  new  high  critical  temperature 
superconducting  materials  that  may  then  be  used  in  beneficial  technological  applications.  Parallel 
to  these  efforts,  it  is  important  to  study  the  macroscopic  behavior  of  superconductors.  Indeed,  once 
new  and  usable  materiads  are  identified,  their  incorporation  into  superconducting  devices  depends 
on  being  able  to  model  these  devices  and  to  numerically  simulate  their  physical  behavior.  Thus, 
scientists,  engineers,  and  mathematicians  who  wish  to  design  superconducting  devices  and/or  who 
wish  to  study  the  physics  of  superconductivity  are  in  need  of  robust  and  efficient  algorithms  and 
codes  for  the  numerical  simulation  of  superconducting  phenomena.  Our  goals  are; 

•  to  develop,  analyze,  and  implement  algorithms  that  are  applicable  to  high-temperature 
superconductors;  and 

•  to  use  these  codes  to  study,  in  coUaboration  with  physicists  and  material  scientists,  super¬ 
conducting  phenomena. 

We  have  already  passed  a  necessary  stepping  stone  towards  meeting  these  goals,  namely  the  devel¬ 
opment  and  implementation  of  methods  for  low- temperature  superconductors.  In  addition,  we  have 
obtained  new  analytical  results  in  this  setting.  A  summary  of  the  results  obtained  are  reported  on 
in  the  papers  listed  in  §III.  Our  work  on  superconductivity  is  receiving  widespread  recognition.  We 
have  been  invited  to  numerous  meetings  and  workshops  on  the  subject,  and  physicists  at  government 
labs,  universities,  and  industry  have  expressed  interest  in  collaborating  with  us  and  using  our  algo¬ 
rithms.  One  collaborative  program  that  has  been  initiated  is  with  the  group  at  Oxford  University 
headed  by  John  Ockendon. 

The  major  accomplishments  resulting  from  our  work  on  low-temperature  superconductivity  is 
that  our  work  represents  the  first  successful  simulation  of  superconductivity  by  a  standard  numerical 
technique  (a  finite  element  method);  our  algorithms,  compared  to  previous  efforts,  are  efficient, 
robust,  and  extendible  to  more  complex  settings,  including  some  involving  high  critical  temperature 
superconductors. 


2.1.  Summary  of  theoretical  results 

We  summarize  the  theoretical  aspects  of  our  preliminary  work;  most  details  can  be  found  in 
the  papers  on  superconductivity  listed  in  §11. 

We  first  focused  on  Ginzburg-Landau  models  for  bounded  two-  or  three-dimensional  regions 
representing  material  samples.  We  reviewed  the  mathematical  formulations  of  some  important 
physical  concepts  such  as  the  fluxoid  c  iantization  emd  several  important  scales  and  parameters. 
Mathematically,  the  electromagnetic  state  of  the  superconductor  corresponds  to  a  minimizer  of  the 
Gibbs  free  energy  which,  in  non-dimensionalized  form,  may  be  formulated  as 


-I-  ih|2  -  2h  •  H  dfi , 


(7) 


where  A,  and  h  =  curl  A  denote  the  non-dimensionalized  complex  order  parameter,  magnetic 
potential,  and  magnetic  field,  respectively;  H  is  the  applied  field  he  constant  /„  is  the  free  energy 
density  of  the  normal  state  in  the  absence  of  a  msgnetic  field,  and  k,  known  as  the  Ginzburg-Landau 
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parameter,  is  a  material  constant  representing  the  ratio  of  penetration  length  to  the  coherence  length. 
For  type-I  superconductors,  we  have  k  <  l/\/2,  while  k  >  l/^/2  for  type-II  superconductors.  A 
formal  study  on  the  transition  between  diflFerent  slates  in  both  the  type-I  and  type-II  superconductors 
has  also  been  carried  out. 

The  Ginzburg-Landau  functional  (7)  has  a  very  important  property,  namely,  that  of  gauge 
invariance.  That  is,  if  for  some  (i>  €  the  standard  Sobolev  space  whose  elements  and  their 

derivatives  of  order  up  to  two  are  square  integrable  on  fl,  we  have 

C  =  and  Q  =  A  4-  Vc2> , 


then  (<,■,  Q)  and  (i/»,  A)  are  said  to  be  gauge  equivalent,  note  that  Q)  =  Q{xi>,  A).  Based  on  this, 
the  existence  of  minimizers  of  the  Ginzburg-Landau  functional  in  the  space  H^{Q)  x  H^{Q)  was 
established  via  standard  variational  arguments.  (For  definitions  of  the  related  spaces,  we  refer  to  the 
papers  listed  in  §III.)  In  fact,  we  showed  that  any  miuimizer  of  the  Ginzburg-Landau  functional  is 
gauge  equivalent  to  a  solution  which  has  a  divergence  free  magnetic  potential  with  vanishing  normal 
component  on  the  boundary.  This  corresponds  to  the  Coulomb  gauge.  The  Ginzburg-Landau 
equations  and  natural  boundary  conditions  in  the  Coulomb  gauge  are  given  by 


2 

—  \p  -r-  =  0 


in  fl , 


(8) 


curlcurl  A  = -  ti'Vii'')  -  Ivl^A ccrl  H  inQ,  (9) 

2k 

div  A  =  0  in  fl  and  A  •  n  =  0  on  F  ,  (10) 

and 

i 

(— Vti)  4- Ai^)  •  n  5=  0  and  curlAxa  =  Hxn  on  F  .  (11) 

K 

Other  {jossible  boundary  conditions  were  also  considered.  Various  mathematical  properties  of  the 
solutions  of  the  Ginzburg-Landau  equations  were  studied.  These  included  the  proof  of  non-existence 
of  local  maxima  for  the  Ginzburg-Landau  functional  and  the  fact  that  the  order  parameter  is  bounded 
by  its  ideal  superconducting  value,  i.e.,  in  non-dimensionaiized  form,  \ii\  <  1  almost  everywhere. 
Simple  analytic  solutions  and  trivial  solutions  were  discussed.  Many  of  the  properties  may  be  used 
to  interpret  related  physical  phenomena  such  as  the  perfect  Meissner  effect  in  the  absence  of  an 
applied  field  and  the  existence  of  a  mixed  state,  i.e.,  the  existence  of  vortices  or  filaments,  in  the 
presence  of  an  applied  field  below  a  critical  value. 

We  also  developed  and  analyzed  finite  element  algorithms  for  approximating  solutions  to  the 
model.  Finite  element  subspaces  are  constructed,  in  a  standard  way,  from  partitions  of  Q  into  finite 
elements;  h  will  denote  some  measure  of  the  size  of  the  finite  elements  in  the  psurtition.  Finite  clement 
approximations  are  then  defined  based  on  a  weak  formulation  of  the  Ginzburg-Landau  equations 
with  the  gauge  div  A  =  0  in  fl  and  A  •  n  =  0  on  T.  We  proved  the  convergence  of  the  finite 
element  solutions  to  a  branch  of  regular  solutions  of  the  nonlinear  Ginzburg-Landau  equations  and 
the  convergence  was  shown  to  be  uniform  for  in  a  compact  interval.  Optimal  error  estimates  were 
also  derived  under  the  usual  regularity  assumptions,  ’however,  we  note  that  the  question  of  the 
regularity  of  the  solutions  of  the  Ginzburg-Landau  equations  in  bounded  domains  has  not  been  fully 
resolved. 
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While  our  analyse*  an4  algorichiiis  were  vaiij  for  hoth  ibe  '.yptj-i  30*4  type- 11  regsiiies,  .i  is  <.4  uar 
moatiy  tor  type-1  supercou4uetors-  It  ji  well-known  that  -a  type-ll  juperoonOucutrs  the  suiutjuiw 
exhibit  much  more  complicated  structures.  The  griti  •nm  tjecessaary  to  jewoh-e  the  %m  struct uf’es 
baaed  on  bounded  domain  modeia,  m  any  computatiots  of  praiticai  ouhty.  wouia  !>e  proioiauveyv 
small.  So,  we  turned  our  attention  to  a  perstMin:  dtfirhury'l..onuuu  mod**  wiucii,  hum  a  imuciu-su 
viewpoint,  is  more  suitable  for  the  aumertcal  Simulation  of  tyj-w*-!!  superfonductofs 

The  periodic  model  is  based  on  trie  uuta>a  of  .Huxoid  -quant. r.at.on  and  -s  ;orMiUi»i,eii  •,*.»  .'.eswriixei 
two-dimensional  thin  rilms.  The  central  asaumptiun  ;n  trie  ierivatioo  of  the  nsotSei  .»  i.ha.i  awas 
from  bounding  surfaces,  certain  physical  variabien.  such  a*  the  ma^nriic  held,  the  cunent  and 
the  density  of  superconducting  charge  earners,  exhibit  penmhc  t>ehav;ur  with  rewpwrt  to  a  tm>- 
dimensional  lattice;  the  accompanying  Sgure  shows  a  typical  ceil  ;a  the  latSi  >  The  aaLin-orthtiguriai 
lattice  is  not  necessarily  aiignesl  with  the  ..xe.irdinate  axes  conauit  tiie  hgure-  A  functioe  /  ,s 
penodic  with  nspect  to  the  lattice  <je<ctori  tj  and  tj  if 

f'x  -  t*;  -  /iX;  for  «  ^  1.  J  and  U'"  ^  12': 


The  cel!  ftp  determined  by  the  lattice  vectors  tj  and  and  the  point  P 

We  were  first  concerned  with  the  formulation  of  the  mcKlei.  A  very  important  iasue  i.s  the 
choice  of  gauge.  A  gauge  choice  is  constrained  by  the  ne<si  to  enforce  she  periodicity  of  several 
physical  variables  and  by  the  fiuxoid  quantization  condition;  it  is  made  more  difficult  by  the  fact 
that  these  constraints  result  in  nonstandard  ■'quasi" -periodic  conditions  for  the  primary  variables 
in  the  Ginzburg-Landau  model,  i.e.,  the  order  parameter  and  magnetic  potential.  There  have  been 
conflicting  discussions  related  to  this  issue  in  the  literature.  W'e  made  a  rigorous  and  detailed  study 
and  were  able  to  define  consistent  gauge  choices.  In  short,  given  the  lattice  cell  depicted  in  the  above 
figure,  we  showed  that  the  vector  magnetic  potential  is  gauge  equivalent  to  a  potential  of  the  form 
Q  -  Auo  where 

A„. !(_';)  m 

and  Q  is  divergence  free,  periodic,  and  is  uniquely  determined  up  to  an  additive  constant  vector.  In 
(13h  B  is  the  average  magnetic  field.  To  maintain  the  periodicity  of  the  desired  physical  variables. 
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the  order  parameter  should  satisfy 

C(x  s- 1*)  =  X  €  r_»,  e  =  i,  2  ,  <14} 

where  s?i(x)  =  —B(x  x  t*)/2,  k  ~  1,2,  and  where  r_j  and  r„2  denote  the  left  and  sx>ttoin  ed^gea 
of  the  lattice  ceU,  respectively.  Clearly,  (14)  implies  the  penutiicity  of  the  magnitude  of  the  order 
parameter;  however,  the  phase  suffers  a  jump  across  the  lattice.  The  tiuxoid  quantiCatjon  condition 
holds  under  this  gauge  choice. 

We  discussed  the  equivalence  of  the  liiSerent  forms  of  the  Giazburg-Landau  functional  m  the 
periodic  setting  and  proved  the  existence  of  the  mimmizer  in  the  appropriate  spaces.  Mimmizcrs 
satisfy  a  system  of  equations  similar  to  iS)  and  (9),  with  suitable  modifications  to  the  boundan' 
conditions.  Then,  we  derived  various  properties  of  the  solutions  of  the  Ginzburg-Landau  equations. 
Some  of  the  results,  such  as  the  boundedness  of  the  order  parameters  and  the  discussion  on  simple 
analytic  solutions,  are  similar  to  those  we  obtained  for  the  bounded  domain  case.  Unlike  that  case 
for  which  the  study  of  the  regularity  of  solutions  was  limited,  we  were  able  to  obtain  extensive 
regularity  results  for  the  periodic  problem. 

We  then  considered  finite  element  approximations  of  solutions  of  the  p^-nodic  model.  Penodic 
models  have  been  used  in  the  past  as  a  setting  for  analyzing  and  approximating  phenomena  ;n  ty-pe- 
II  superconductors  using  for  the  most  part,  some  type  of  senes  solution  or  a  Monte  Carlo/simulated 
annealing  approach.  .-Vs  in  other  settings,  finite  element  methods  can  be  very  competitive  for  nu¬ 
merical  simulation  purposes. 

Basically,  we  employ  the  standard  Galerkin  finite  element  approximation  to  solve  the  nonlinear 
Ginzburg-Landau  eqimtions.  The  periodicity  of  the  physical  variables  and  boundary  conditions  such 
as  (14)  icr  the  primary  variables  a*e  the  non-standard  features  of  our  scheme.  In  general,  functions 
in  the  finite  element  space  only  satisfy  (14)  at  the  interpolating  nodes  on  the  boundary,  i.e.. 

+t*)  =  Vx;  e  r_i:.  fc  =  1.2.  (15) 

where  Xj  is  any  boundary  vertex  of  a  triangle  for  piecewise  linear  and  quadratic  elements  and  it  may 
also  be  the  midpoint  of  any  triangle  edge  on  the  boundary  for  quadratic  elements.  The  periodicity 
of  finite  element  solutions  with  respect  to  the  lattice  is  also  defined  only  at  interpolating  nodes. 
Thus,  the  finite  element  spaces  are  not  the  subspaces  of  the  exact  solution  spaces.  The  study  of  the 
approximations  centered  around  this  issue. 

The  key  to  the  error  analysis  is  an  estimate  of  a  boundary  integral  term  which  would  disappear  if 
the  finite  element  spaces  were  subspaces  of  the  underlying  solution  space  for  the  continuous  problem. 
We  showed  that  this  integral,  although  not  zero  in  general,  gives  a  higher  order  error  term  in  the  final 
estimate  and  therefore,  the  optimal  convergence  rates  are  retained.  The  idea  may  well  be  generalized 
to  higher  order  element  cases  with  suitable  choices  of  interpolation  procedures.  Thus,  we  were  able 
to  obtain  optimal  error  estimates  for  approximations  to  the  solution  of  the  full  nonlinear,  periodic, 
Ginzburg-Landau  equations. 

2.2.  Summary  of  computational  results 

We  have  developed  a  finite  element  code  for  the  periodic  Ginzburg-Landau  model  discussed 
above.  The  code  uses  piecewise  quadratic  finite  element  functions  based  on  a  triangulation  of  a 
lattice  cell.  Tlie  periodicity  and  “quasi” -periodicity  conditionr  eure  implemented  as  described  by 
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(15),  i.e.,  at  the  aodes  of  the  finite  ekmenta.  The  moet  interesting  periodicity  structure  is  that  of 
an  equilateral  triangular  Lattice  having  one  fiuxoid  associated  with  each  lattice  ceU.  Indeed,  it  is 
well-known  that  such  a  lattice  yields  the  smallest  value  for  the  Gibbs  free  energy.  Having  chosen 
input  data  corresponding  to  an  equilateral  triangular  lattice  having  one  fluxoid  associated  with 
each  lattice  ceil,  the  only  remaining  inputs  to  be  chosen  are  the  Gin^burg-Landau  parameter  a.  the 
average  magnetic  field  B,  and  the  number  and  spacing  of  the  grid  lines. 

The  discretized  equations  are  a  nonlinear  system  of  algebraic  equations.  These  are  solved  by 
a  continuation  method  coupled  with  Newton’s  method.  The  code  is  configured  so  that  one  chooses 
a  fixed  value  for  k,  and  then  varies  5.  For  each  pair  (a,  5),  Newton’s  method  is  used  to  solve  the 
nonlinear  equations.  The  initial  guess  for  Newton’s  method  is  determined  by  continuing  from  the 
solution  determined  for  a  previous  value  of  5  and  the  same  value  of  k.  The  particular  continuation 
method  used  is  a  tangent  line  approximation  to  the  solution  at  the  previous  value  of  B.  We  start 
with  a  value  of  B  close  to  the  upper  limit  «  for  which  Newton’s  method  seems  to  have  a  large 
attraction  bail;  we  then  continue  by  successively  reducing  the  value  of  B  towards  its  lower  limit  0. 

The  graph  of  the  computed  approximation  of  the  magnetization  -AirSi  vs.  the  applied  field  He 
and  the  level  curves  of  the  density  of  superconducting  electrons  .V,  are  given  in  the  next  two  figures, 
respectively.  These  results  were  obtained  on  a  Macintosh  II  using  3  uniformly  spaced  intervals  in 
each  of  the  lattice  directions.  For  comparison  purposes,  we  also  provide,  in  the  first  figure,  the 
corresponding  graph  fir  a  Monte  Carlo/simulated  annealing  approximation.  For  the  second  figure, 
the  solution  in  only  a  single  lattice  cell  was  computed;  this  solution  was  extended,  using  periodicity 
or  “quasi” -periodicity  relations  to  obtain  the  solution  outside  the  computational  cell. 


Comparison  of  computational  results  for  the  magnetization  (-4n'M)  vs.  external 
field  (^fe)  for  «  =  5. 


Other  pa$t  actixhttes 
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Level  curves  of  Ns,  the  density  of  superconducting  charge  carriers,  for  «  =  5  and 
S  =  10/3. 


II  -  SUMMARY  OF  OTHER  RESEARCH  PROJECTS 


Under  the  sponsorship  of  AFOSR  Grant  Number  AFOSR-90-0179,  we  have  also  engaged  in  a 
number  of  other  research  projects.  Here  we  give  brief  descriptions  of  some  of  our  accomplishments. 
Details  may  be  found  in  the  papers  listed  in  §III. 

•  Parallel  algorithms  for  flow  calcidations  based  on  the  velocity-vorlicity  formulation.  We 
carried  out  a  careful  study  of  the  accuracy  of  the  approximations  of  the  Navier-Stokes 
through  the  use  of  the  velocity-vorticity  formulation.  We  have  shonw  that  most  of  the 
existing  boundary  treatments  lead  to  bad  vorticity  approximations,  and  we  have  also  de¬ 
veloped  improved  treatments.  We  have  demonstrated  how  this  formulation  can  be  used  to 
advantage  in  a  parallel  computing  envirnroent.  Specifically,  for  methods  of  solution  of  the 
Navier-Stokes  equations  that  require  the  solution  of  a  series  of  Stokes  problems,  we  have 
shown  how  the  latter  may  be  solved  as  a  sequence  of  uncoupled  Poisson  problems.  Thus 
we  have  coarse-grain  parallelism;  fine-grain  parallelism  may  also  be  attained  within  each 
of  the  Poisson  solves. 

•  Analysis,  application,  and  computation  of  centroidal  Voronoi  tesselations.  Voronoi  tesse- 
lations,  and  their  dual  Delauny  tesselations  of  great  use  in  a  variety  of  seetings,  including 
numerical  analysis,  data  compression,  etc.  We  examine  the  special  type  of  Voronoi  tessela- 
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tion  wherein  the  points  from  which  the  Voronoi  regions  axe  denned  are  also  the  centroids 
of  thtae  regions.  We  have  shown  that  this  type  of  tesaelation  possess  a  certain  optimality 
property  that  makes  them  useful  in  applications.  We  have  also  analyzed  algorithms  for  the 
determination  of  such  tesselations. 

*  The  treatment  of  irhomogeneons  essenttai  boundary  conditions  and  the  accurate  computa¬ 
tion  of  stresses  along  boundaries  in  the  Jintte  element  method.  The  ■‘basic"  hxute  element 
method  applies  to  homogeneous  Dirichlet  boundary  oonditiotis.  A  variety  of  treatments 
of  inhomogeneous  boundary  conditions  has  been  proposed  in  the  Literature.  In  our  work, 
we  studied  the  moat  popular  treatment  in  engineering  practice,  namely  approximating  the 
inhomogeous  data  by  a  boundary  interpolant.  We  also  studied  Lagrange  multiplier  meth¬ 
ods,  and,  in  particular,  developed  a  claaa  of  specific  methods  -hat  are  optiraaily  accurate 
and  for  which  the  Lagrange  multiplier  computation  uncouples.  We  made  a  comprehensive, 
rigorous  study  of  the  accuracy  of  these  various  treatments.  As  a  side  benefit,  we  devised  an 
algorithm  for  the  accurate  calculation  of  boundary  stresses  in  fluid  flow  calculations.  We 
also  analyzed  the  accuracy  of  these  stress  approximations. 


Ill  -  SUMMARY  OF  OTHER  ACTIVITIES  UNDER  GRANT  AFOSR-90-0179 


Book  editing  project 

•  Editing  a  book  (with  RA.  Nicolaides)  to  be  published  by  Cambridge  University  Press  on 
recent  trends  and  advances  in  incompressible  flow  calculations. 

Personnel  supported  by  the  grant 

•  Max  Gunzburger  (Principal  investigator)  -  Summer  salary. 

•  Jerome  Eastham  (Visiting  Assistant  Professor)  -  Partial  academic  year  salary. 

•  Mark  Mundt  (Ph.D.  student)  -  Academic  year  salary. 

•  Pavel  Bochev  (Ph.D.  student)  -  Summer  salary. 

•  H.  Lee  (Ph.D.  student)  -  Summer  salary. 

•  John  Burkhardt  (Ph.D.  student)  -  Academic  year  and  summer  salary. 

Conferences,  workshops,  and  special  sessions  organised  on  grant  related  research 

•  Special  Session  on  Control  of  Partial  Differential  Equations,  World  Congress  of  Nonlinear 
Analysts,  Tampa,  August,  1992,  (with  J.  Burns  and  T.  Herdman). 

•  Workshop  on  Flow  Control,  IMA,  Minneapolis,  November,  1992. 

•  Special  Session  on  Flow  Control,  IEEE  Conference  on  Decision  and  Control,  Tuscon,  De¬ 
cember,  1992,  (with  K.  Ito). 

•  AMS-SIAM-IMS  Summer  Research  Conference  on  Superconductivity,  site  to  be  announced, 
July,  1993,  (with  J.  Ockendon). 

Invited  talks  in  1991-1992  on  grant  related  research 

•  Internationed  Conference  on  Differential  Equation,  Edinburg,  TX,  May,  1991. 

•  Seventh  International  Conference  on  Numerical  Methods  in  Laminar  amd  Turbulent  Flow, 
Stanford,  July,  1991. 

•  Second  Soviet-North  American  Workshop  on  Computational  Fluid  D)mamics,  Montreal, 
September,  1991. 

•  Fourth  International  Symposium  on  Computational  Fluid  Dynamics,  Davis,  CA,  Septem¬ 
ber,  1991. 

•  Workshop  on  Superconductivity,  Argonne,  IL,  January,  1992. 

•  The  Mathematics  of  Superconductivity,  Edinburgh,  Scotland,  January,  1992. 

•  Joint  Meeting  of  the  American  and  London  Mathematical  Societies,  Cambridge,  England, 
June,  1992. 

•  World  Congress  of  Nonlinear  Analysts,  Tampa,  August,  1992. 


Other  pa»t  acttvtties 
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•  Computation  and  Control  III,  Bozeman,  MT,  .August,  1992. 

•  Sixth  International  Conference  on  Boundary  and  Interior  Layers,  Summit  Country,  CO. 
August,  1992. 

•  Conference  on  Computational  Methods  for  the  Material  Sciences,  Pittsburgh.  September, 
1992. 

•  American  Mathematical  Society  Meeting,  Dayton,  October,  1992. 

•  IEEE  Conference  on  Decision  and  Control,  Tuscon,  December,  1992 

Piiblications  prepared  reporting  on  grant  related  researth  (]  990-1 992) 

1.  Finite  element  approximations  of  a  Ladyzhenskaya  model  for  stationary  incomprtsKibie 
viscous  flow;  SIAM  J.  Numer.  Anai.  2T,  1990,  1-19;  with  Q.  Du. 

2.  Elxperiences  with  computational  methods  for  the  velocity-vorticity  formulation  of  incom¬ 
pressible  viscous  flow;  Computational  Methods  in  Viscous  .Aerodynamics,  Elsevier.  1990, 
231-271;  with  M.  Mundt  and  J.  Peterson. 

3.  A  subdomain  Galerkin/least  squares  method  for  first  order  elliptic  systems  in  the  plane; 
SIAM  J.  Numer.  Anal.  27,  1990,  1197-1211;  with  C.  Chang. 

4.  A  numerical  method  for  drag  minimization  via  the  suction  and  injection  of  mass  through 
the  boundary;  Stabilization  of  Flexible  Structures,  Springer.  1990,  312-321;  with  L,  Hou 
and  T.  Svobodny. 

5.  Finite  element  approximations  of  an  optimal  control  problem  associated  with  the  scalar 
Ginzburg-Landau  equation;  Comput.  Math.  Appi  21,  1991,  123-131;  with  L.  Hou  and  T. 
Svobodny. 

6.  .Analysis  of  a  Ladyzhenskaya  model  for  incompressible  viscous  flow;  J.  Math.  Anai  Appi 
155,  1991,  21-45;  with  Q.  Du. 

7.  Analysis,  approximation,  and  computation  of  control  problems  for  incompressibe  flows; 
Turbulence  Structure  and  Controi  Ohio  State,  1991,  85-88. 

8.  Existence,  uniqueness,  and  finite  element  approximation  of  solutions  of  the  equations  of 
stationary,  incompressible  magnetohydrodynamics;  Math.  Comp.  56,  1991,  523-563;  with 
A.  Meir  and  J.  Peterson. 

9.  Control  of  temperature  distributions  along  boundaries  of  engine  components;  Numerical 
Methods  in  Laminar  and  Turbulent  Flow  VII,  Pineridge,  1991,  765-773;  with  L.  Hou  and 
T.  Svobodny. 

10.  Vorticity  constraints  in  velocity-vorticity  formulations  of  steady,  viscous,  incompressible 
flow;  Numerical  Methods  in  Laminar  and  Turbulent  Flow  VII,  Pineridge,  1991,  774-781; 
with  Q.  Du  and  A.  Meir. 

11.  Approximation  of  boundary  control  and  optimization  problems  for  fluid  flows;  lift  Interna¬ 
tional  Symposium  on  Computational  Fluid  Dynamics,  U.  California,  Davis,  1991,  455-460; 
with  L.  Hou  and  T.  Svobodny. 

12.  Analysis  and  finite  element  approximations  of  optimal  control  problems  for  the  stationary 
Navier-Stokes  equations  with  distributed  and  Neumann  controls;  Math.  Comp.  57  1991, 
123-151;  with  L.  Hou  and  T.  Svobodny. 

13.  Predictor  and  stepiength  selection  in  continuation  methods  for  the  Navier-Stokes  equations; 
Comput.  Math.  Appi  22,  1991,  73-81;  with  J.  Peterson. 

14.  Analysis  and  finite  element  approximations  of  optimal  controi  problems  for  the  stationary 
Navier-Stokes  equations  with  Dirichlet  controls;  Math.  Model  Numer.  Ancd.  mbf  25, 
1991,  711-748;  with  L.  Hou  and  T.  Svobodny. 

15.  Boundary  velocity  control  of  incompressible  flow  with  an  application  to  viscous  drag  re¬ 
duction;  SIAM  J.  Control  Optim.  30  1992,  167-181;  with  L.  Hou  and  T.  Svobodny. 

16.  Analysis  and  approximation  of  the  Ginzburg-Landau  model  of  superconductivity;  SIAM 
Review  34  1992,  54-81;  with  Q.  Du  and  J.  Peterson. 

17.  Numerical  solution  of  the  compressible  boundary  layer  equations  using  the  finite  element 
method;  AIAA  Paper  AIAA-92-0666,  AlAA,  Washington,  1992;  with  E.  Hjtopouios  and 
J.  Schetz. 


Other  past  acJ>tnJ»ej 


18.  IVeating  inhomc^eneous  essential  boundary  conditions  in  baile  element  methods  and  the 
calculation  of  boundary  stresses;  SIAM  J.  Numer.  Anal  29  1992,  390-424;  with  L.  Hou. 

19.  Solving  the  Ginzburg- Landau  equations  by  finite  element  methods;  Phys.  Rev.  B-.  46  1992, 
9027-9034;  with  Q.  Du  and  J.  Peterson. 

20.  On  the  Ginzburg-Landau  equations  of  superconductivity;  Partial  Differenttal  Equations, 
Longman,  1992,  58-62;  with  Q,  Du  and  J.  Peterson. 

21.  Heating  and  cooling  control  of  temperature  distributions  along  boundaries  of  flow  domains; 
to  appear  in  J.  Math.  Systems  Esttm.  Control  with  L.  Hou  and  T.  Svobodny. 

22.  Optimal  control  and  optimization  of  viscous,  incompressible  flows;  to  appear  in  Incompress¬ 
ible  Computational  Fluid  Dynamics:  Trends  and  Advances;  with  L.  Hou  and  T.  Svobodny. 

23.  Boundary  control  of  incompressible  flows;  to  appear  in  Advances  in  Computational  Fluid 
Dynamics,  with  L.  Hou  and  T.  Svobodny. 

24.  Modeling  and  analysis  of  a  periodic  Ginzburg-Landau  model  for  type-II  superconductors; 
to  appear;  with  Q.  Du  and  J.  Peterson. 

25.  Optimal  control  problems  for  a  class  of  nonlinear  equations  with  an  application  to  the 
control  of  fluids;  to  appear;  with  L.  Hou  and  T.  Svobodny. 

26.  Finite  element  approximation  of  a  periodic  Ginzburg-Landau  model  for  type-II  supercon¬ 
ductors;  to  appear;  with  Q.  Du  and  J.  Peterson. 

27.  The  approximation  of  boundary  control  problems  for  fluid  flows  with  an  application  to 
control  by  heating  and  cooling;  to  appear;  with  L.  Hou  and  T.  Svobodny. 

28.  Analysis,  applications,  and  computation  of  centroidal  Voronoi  tesselations;  in  preparation; 
with  Q.  Du,  V.  Faber,  and  C.  Scovill. 


